Could 2S 0114+650 be a magnetar? 
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ABSTRACT 

We investigate the spin evolution of the binary X-ray pulsar 2S 0114+650, which possesses 
the slowest known spin period of ~ 2.7 hours. We argue that, to interpret such long spin period, 
the magnetic field strength of this pulsar must be initially £ 10 14 G, that is, it was born as a 
magnetar. Since the pulsar currently has a normal magnetic field ~ 10 12 G, our results present 
support for magnetic field decay predicted by the magnetar model. 

Subject headings: binaries: close - pulsars: individual: 2S 0114+650 - stars: neutron - X-ray: 
stars 



Introduction 



Neutron stars are thought to be born as rapidly rotating (~ 10 ms) radio pulsars created during a 
type Il/Ib supernova explosion involving a massive star. The dipolar magnetic fields of radio pulsars, as 
inferred from their observed spin-down rates, range from 10 8 G to 3 x 10 13 G (Taylor, Manchester, & Lyne 
1993). However, it has been proposed that there may exist "magnetars" - neutron stars with magnetic 
field strengths in excess of ~ 10 14 G (Thompson & Duncan 1992). These objects have been invoked to 



model soft gamma-ray repeaters (SGRs) and the 6 — 12 s anomalous X-ray pulsars (AXPs) ( Thompson 
Duncan 199& see also |Kouvcliotou ct al. 19931 , 1994; [van Paradijs, Taam, fc van den Heuvel 1995^ |Corbet 



ct al. 1995 ; Vasisht, Frail, fc Kulkarni 1995 ; Vasisht fc Gotthelf 1997 ), though unambiguous evidence for 
the existence of magnetars comes from recent observations of the AXP-like object IE 1841—045 (Vasisht & 



Gotthelf 1997), SGRs 1806-20 (Kouvcliotou ct al. 1998a) and 1900+14 (Kouveliotou et al. 19981) 



According to Thompson & Duncan (1992, 1996), magnetars are neutron stars born with millisecond 
periods that generate magnetic fields above 10 14 G by dynamo action due to convective turbulence, 
magnetic field decay powers the X-ray and particle emission of magnetars. It is conceivable that when the 
magnetic fields in magnetars have decayed to, say, 10 12 G, and accretion occurs, they are not distinguished 
from X-ray pulsars born with normal (10 12 — 10 13 G) magnetic fields, except that they may have relatively 
longer spin periods. One may then expect to find magnetar descendants among long-period binary X-ray 
pulsars, because they are much more luminous in X-rays than isolated objects accreting from interstellar 
medium. Here we present arguments indicating that 2S 0114+650, the X-ray pulsar with the slowest known 
spin, may have a magnetar evolutionary history. 



2. The slowest X-ray pulsar 2S 0114+650 



The X-ray source 2S 0114+650 was discovered in 1977 by the SAS 3 galactic survey (Dower fc Kelly 



1977). Its optical counterpart, LS 1+65 010, was recently identified as a supergiant of spectral type Bl la 
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(Rcig et al. 1996). Thus 2S 0114+650 belongs to the class of high-mass X-ray binaries (HMXBs), systems 
in which a compact star - generally a neutron star - accompanies a high-mass donor star (cf. Bhattacharya 
fe van den Heuvel 1991 for a review). With a distance of 7.2 kpc derived from this spectral classification, 
the X-ray luminosity is a few 10 35 — 10 36 ergs _1 . An orbital period of 11.59 days was reported from optical 
radial velocity measurements by Crampton, Hutchings, & Cowley (1985). There has been some weak 
evidence of a pulsation period of 850 — 895 s ( Yamauchi et al. 1990] ; Koenigsberger et al. 1983| ). In contrast, 
Finlcy, Belloni, & Cassinelli (1992) have discovered periodic X-ray outbursts with a 2.78 hour period. The 
same period was confirmed by ROSAT observations ( Finley, Taylor, fe Belloni 1994 ). Recent Rossi X-ray 
Timing Explorer (RXTE) observations show presence of modulation on 11.59 day orbital period as well as 
2.7 hour pulse period ( Corbet, Finley, fc Peele 1998 ). If this pulse period indeed represents the rotation 
period of the neutron star, 2S 0114+650 would be by far the slowest known X-ray pulsar. 



3. The spin evolution in 2S 0114+650 



How has 2S 0114+650 been spun down to the long period (P ~ 10 4 s) if it was formed with much 
shorter period (~ 0.01 — Is, say)? The neutron star's spin evolution is divided in three phases (see Davies 



fc Pringlc 1981; Bhattacharya & van den Heuvel 1991). In the first radio pulsar phase, the star is an 



active radio pulsar and spins down by magnetic dipole radiation. This phase ends when the ram pressure 
of the ambient material overcomes the pulsar's wind pressure at the gravitational radius (Rq = 2GM/V 2 , 
where G is the gravitation constant, M is the neutron star mass, and V is the relative velocity between 
the neutron star and the ambient material). In the second propeller phase material enters the corotating 
magnetosphere and is stopped at Ra, the Alfven radius, where the energy density in the accretion flow 
balances the local magnetic pressure. Further penetration cannot occur owing to the centrifugal barrier; 
that is, Ra > R c , where R c = [GM(P/2n) 2 ] 1 / 3 is the corotation radius. The star expels the material once 
it spins up the material to the local escape velocity at ~ Ra (Illarionov & Sunyaev 1975). Propeller action 
continues until Ra — R c , when the centrifugal barrier is removed and the spin period reaches its equilibrium 



value (Bhattacharya fc van den Heuvel 1991) 



P oq * (20s)B% 7 Mj /7 Rf /7 M 1 
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(1) 



where B = 10 Byi G is the neutron star's dipolar magnetic field strength, M = 10 Misgs the 
mass accretion rate, R — 10 6 i?6 cm the radius, and Mi. 4 = M/1.4M (Throughout this Letter, we take 
-Mi. 4 = Rq = 1). In the following accretion phase the star becomes an X-ray pulsar, and its spin evolution 
is determined by the net torque exerted on the star by the accretion flow. 

The X-ray emission observed in 2S 0114+650 is most likely to be powered by accretion onto the neutron 
star via a stellar wind from the companion. With a mass-loss rate of a few 10 -6 M Q yr _1 and a wind velocity 
of ~ 10 kms , the derived X-ray luminosity from a simple wind- fed model is in accordance with the 

1996| ). A stable accretion disk is unlikely to form around the neutron star, 

Early two-dimensional 



200 km s" 1 ) wind velocity (Wang 1981 



mean detected one ( Reig et al. 
which would require an extremely low (<- 
numerical studies of Bondi-Hoyle accretion flow (e.g., Matsuda, Inouc, fc Swada 1987 : Fryxcll fc Taam 1988] , 
1989) demonstrated that temporary accretion disks with alternating sense of rotation possibly form in a 



wind accreting system. More recent high resolution three-dimensional numerical investigations (e.g., Ruffcrt 
1992) . 1997) found the so-called wind "flip-flop" instability with the timescale of the order of hours. This is 
consistent with the torque fluctuations in wind-accreting X-ray pulsars (Nagase 1989), suggesting that the 
long-term averaged angular momentum transferred to the neutron star by the accreted wind material is 
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very small. We therefore conclude that the spin period of 2S 0114+650 has not been considerably changed 
during the present accretion phase. This means that the long spin period of 2S 0114+650 was actually 
attained in an earlier evolutionary phase before the companion star became a super-giant - it was spun down 



by the propeller mechanism (Qlarionov & Sunyacv 1975) when the companion star was on main-sequence 
and had a weaker wind. The magnitude of the equilibrium period, as seen in equation (1), is determined by 
the magnetic field strength and mass accretion rate of the neutron star during this phase. 

The magnetic field strength in 2S 0114+650 has not been measured, since no cyclotron features have 
been seen in its X-ray spectrum. But there is indirect evidence indicating that its magnetic field strength 
is similar to those in other X-ray pulsars: Its spectrum shows the typical shape of the usual X-ray pulsars 
having a power-law with an exponential high-energy cutoff at ~ 7 keV (Yamauchi et al. 1990) or > 15 



keV (Koenigsbcrgcr ct al. 1983). The iron emission line at about 6.4 keV, which is common among X-ray 



pulsars, was also discovered in 2S 0114+650 (Yamauchi et al. 1990; Apparao, Bisht, & Singh 1991) 



Observations of cyclotron lines in X-ray pulsars imply surface fields of about (0.5 — 5) x 10 12 G, and it has 
been suggested that the cutoff seen in the power-law spectra of X-ray pulsars is related to the magnetic 
field strength of the neutron star ( Makishima fc Mihara 1992] ) ^|. This would imply a field of ~ 10 12 G for 
2S 0114+650, consistent with those for typical X-ray pulsars which show cutoff energies of 10 — 20 keV. 



The mass of LS 1+65 010 was estimated to be 16(±5)M© from evolutionary models ( Reig et al. 1996 ) 
For a typical mass-loss rate of ~ 10~ 8 — 10 _7 M Q yr _1 from a 16 M & main-sequence star and a wind 



velocity of ~ 10 3 kms , the accretion rate of the neutron star is M 



10 



13 



lO^gS" 1 . 



Inserting the 



values of B and M into equation (1), we find P cq ~ 50 — 140 s, nearly two orders of magnitude shorter than 
the observed period. (Even if B is enhanced to 10 13 G, P eq never exceeds ~ 1000 s.) This implies that the 
extremely long period of 2S 0114+650 cannot be reached via the propeller mechanism if the pulsar has 
possessed a constant magnetic field of ~ 10 12 — 10 13 G. One may argue that the accretion rate during the 
propeller phase could be much lower than that adopted here, because of a lower rate of mass loss from the 
companion star or a higher wind velocity. For example, if M ranges from ~ 10 9 gs _1 (for B = 10 12 G) to 
~ 10 11 gs _1 (for B = 10 13 G), the value of P eq can be indeed raised to ~ 10 4 s. However, this would lead 



to a spin-down time during the radio pulsar phase (Davies & Pringlc 1981) 

-1/2 



t s ~ (2.5 x 10 10 yr)B^M 9 ~'~V s 
~ (2.5 x 10 8 yi)B^M^ /2 V 8 -1 

(where V — 10 8 Vs cms -1 ), which is much longer than the lifetime the companion star spends on 
main-sequence (generally 10 7 yr). 



(2) 



There exists another possibility that the neutron star was born rotating slowly (P ~ 1 s), so that it 
went directly to the propeller phase. Propeller spin-down to P cq takes ( Wang fc Robertson 1985| ) 

~ 1/2 Af~ 3/4 p ~ 3/4 



(1.5 x 10 iu yr)B 



9 Pi 



~ (1.5 x 10 8 yr)B t 3 7 Af n ' P x 

(where Pi = P/l s), which is still much longer than 10 7 yr. 

We are eventually led to the conclusion that 2S 0114+650 must initially have had a magnetic field 
much stronger than its present value, that is, it was born as a magnetar. Distinguished from neutron 



(3) 



3 This relation is not quite accurate (see 
estimate of the magnetic fields. 



Reynolds, Parmar, & White 199J), but it may provide an order of magnitude 
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stars with normal magnetic field strengths, magnetars are radio-quiet, and the decaying magnetic fields 
power the X-ray and particle emission. For a magnetar in a binary system, during the early evolutionary 
phase, the pressure exerted by particle emission exceeds the stellar wind ram pressure at ~ Rq, preventing 
accretion onto the neutron star ( Thompson fc Duncan 1996| ). After t ~ 10 4 — 10 5 years, the star's spin 



period increases to 

P~{lOs)tl /2 B 15 R%Mil /2 (4) 

by magnetic dipole radiation (where t± — t/10 4 yr), the particle luminosity decays, and the wind material 
from the companion star begins to interact with the magnetosphere. Again adopt a mass accretion rate of 
- 10 13 - 10 14 gs _1 , the propeller effect, for a magnetic field of - (1 - 4) x 10 14 G, can comfortably spin 
down the neutron star from ~ 10 s to ~ 10 4 s on a timescale ^ 10 5 years (cf. equations [1] and [3]), i.e., 
before the magnetic field decays significantly (see below). After this long equilibrium period is reached, the 
period of the neutron star would remain close to it, due to inefficient angular momentum transfer during 
the subsequent wind-accretion process. 

Field decay in 2S 0114+650 should be slow enough to garantee the final, long spin period, and fast 
enough to allow a considerable reduction in the field within the lifetime of the system. There exist several 
mechanisms for field decay in non-accreting neutron stars: Ohmic decay, ambipolar diffusion and Hall drift 
(e.g., Goldreich & Reisenegger 1992). Ohmic decay dominates in weakly magnetized (B ^ 10 11 G) neutron 



stars, fields in intermediate strength (B ~ 10 12 — 10 13 G) decay via Hall drift, and intense fields (B £ 10 14 
G) are mainly affected by ambipolar diffusion. For typical values of the characteristic length scale (10 5 
cm) of the flux loops through the outer core and the core temperature (10 8 K), the field-decay timescale 
for ambipolar diffusion through the solenoidal mode in a magnetar is around 3 x 10 5 yr if B ~ 10 14 G 



(Thompson & Duncan 1996). Compared to the age (a few 10 6 yr) of the optical companion of 2S 0114+650, 
this implies that the magnetar has had enough time for its field to decay to its current value (a few 10 12 
G). Similar conclusion can also been found from the detailed calculations by Heyl & Kulkarni (1998). 



Discussion 



SGRs, which are transient gamma-ray sources that undergo repeated outbursts, have been suggested 
to be the prototypes of magnetars. There are four known SGRs, with two (SGRs 1806—20 and 0525—66) 
associated with young supernova remnants ( Kulkarni fc Frail 1993| ; Vasisht ct al. 1994 ). Recently, 
pulsations at a period of 7.47 s and 5.16 s with a spin-down rate of 2.6 x 10~ 3 syr -1 and 3.5 x 10~ 3 syr -1 
were discovered in the persistent X-ray flux of SGRs 1806—20 (Kouveliotou ct al. 1998a) and 1900+14 
(Kouveliotou et al. 1998b), respectively; the magnetic field strengths of (2 — 8) x 10 14 G can be derived if 
the spin-down is due to magnetic dipole radiation. AXPs are another type of high-energy sources that have 
recently joined this group of highly magnetized neutron stars. They are a group of about eight pulsating 
X-ray sources with periods around 6 — 12 s (cf. Stella, Israel, fc Mereghctti 1998 for a review), characterized 

10 36 ergs -1 , and very soft 



by steady spin-down, relatively low and constant X-ray luminosities of ~ 10 35 
X-ray spectra. So far no optical counterpart has been detected. Nearly half of them are located at the 
centers of supernova remnants, suggesting that they are relatively young ( 10 5 years). Dipole magnetic 
fields of 10 14 — 10 15 G can also be derived from the measured spin-down rates, if they are spinning down 
due to dipole radiation torques. The observed X-ray luminosities, spin periods and spin-down rates in SGRs 
and AXPs follow naturally from the magnetar model ( Thompson fc Duncan 1996| ). 



A crucial test of the magnetar model should include the observational evidence of the hypothesized 
magnetic field decay in magnetars. This evidence is most likely to be found in slowly-rotating, binary (such 



- 5 - 



as 2S 0114+650) and isolated (such as RX J0720.4-3125; Haberl et al. 1997; Heyl fc Hcrnquist 199S 



Kulkarni fc van Kerkwijk 1998j ) X-ray pulsars. Figure 1 compares the magnetar candidates (SGRs 1806—20 
and 1900+14 in filled rectangles, and AXPs IE 2259+586, 4U 0142+61, IE 1048.1-5937 and IE 1841-05 
in filled stars) with their possible X-ray pulsar descendants (2S 0114+650 in open triangle and other pulsars 
in open stars) in the magnetic field versus spin period diagram. A schematic view of the magnetar field 
and spin evolution of magnetars is clearly seen. The magnetic fields of SGRs and AXPs are derived from 



their observed spin-down rates (Kouveliotou et al. 1998a, 1998b; Stella, Israel, & Mereghetti 1998 and 



references therein) under the assumption that the spin-down is caused by magnetic dipole radiation. For 
2S 0114+650 we take a canonical magnetic field of 3 x 10 12 G with the possible range of 10 12 — 10 13 G. The 
dashed lines represent the relation between the magnetic field strength and the equilibrium spin period with 
various mass accretion rates (10 18 , 10 15 and 10 12 gs _1 ), from which an initial magnetic field of a pulsar can 
be constrained by the present spin period, given a reasonable estimate of the mass accretion rate during 
the propeller phase. The detailed evolutionary track of a magnetar depends on the particle luminosity, the 
properties of the stellar wind from the companion, the magnetic field and its decay timescale. 
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Fig. 1. — The magnetic field versus spin period diagram for the candidates of magnetars (SGRs in filled 
rectangles and AXPs in filled stars) and the slow X-ray pulsars (in open stars except 2S 0114+650 in open 
triangle) with known magnetic fields. The dashed lines denote the relation between the magnetic field and 
the equilibrium spin period with various mass accretion rates. 



In Fig. 1 we have also plotted the slow (spin periods ^ 100 s) X-ray pulsars (A 0535+26, Vela X— 1, 
GX 1+4, 4U 1907+09, 4U 1538-52 and GX 301-2) with known magnetic field strengths as possible 
descendants of magnetars. Their magnetic field strengths are generally estimated from the observed 



cyclotron line features in X-ray spectra (Mihara & Makishima 1998 and references therein). For GX 1+4 



it is determined from the observational evidence of the propeller effect (Cui 1997). These sources may 
have had a similar evolutionary history as 2S 0114+650, but for them the evidence is not as unequivocal 
as for 2S 0114+650. Their slow spins could be accounted for in terms of the propeller effects with normal 



magnetic field strengths (e.g., Illarionov fc Sunyaev 1975 ; Waters fc van Kerkwijk 198E , see also the dashed 
lines). However, a magnetar model presents an alternative explanation that can not be ruled out. Actually, 
high magnetic field strengths ( £ 10 13 G) have been measured in A 0535+26 and GX 1+4. Statistically, 



if the birth rate of magnetars is of the order of one per millennium (Kouveliotou et al. 1998a), and the 
X-ray lifetime of HMXBs lasts up to ~ 10 5 yr, there may exist in the Galaxy ~ 10 binary X-ray pulsars 
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originating from magnetars. 
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